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S Abstract 

D Starting from an initial classical four dimensional fiat background of the heterotic or type II 
r^, superstrings, we are able to determine at the string one-loop level the quantum corrections 
to the effective potential due to the spontaneous breaking of supcrsymmctry by "geometrical 
^ fluxes". Furthermore, considering a gas of strings at finite temperature, the full "effective 
thermal potential" is determined, giving rise to an effective non-trivial pressure. The back- 
^ reaction of the quantum and thermal corrections to the space-time metric as well as to the 
^ moduli fields induces a cosmological evolution that depends on the early time initial condi- 
tions and the number of spontaneously broken supersymmetries. We show that for a whole 
^ set of initial conditions, the cosmological solutions converge at late times to two qualita- 
O tively different trajectories: They are either attracted to (i) a. thermal evolution similar to 
^ a radiation dominated cosmology, implemented by a coherent motion of some moduli fields, 
or to (a) di, "Big Crunch" non-thermal cosmological evolution dominated by the non-thermal 
part of the effective potential or the moduli kinetic energy. During the attraction to the 
radiation-like era, periods of accelerated cosmology can occur. However, they do not give 
rise to enough inflation (e-fold ~ 0.2) for the models we consider, where J\f > 2 supersym- 
metry is spontaneously broken to Af — 0. 
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1 Introduction 



A sensible theoretical description of cosmology has certainly to take into account both the 
quantum and thermal effects. It is naturally the case for an early phase of the history, 
close to the Planck time, where quantum fluctuations are of same order of magnitude as the 
whole size of the universe and where the temperature is extremely high. It has to be the 
case as well at very late times, where the description of the cosmology should involve the 
electroweak phase transition and provide the missing link to the quantum particle physics 
of the standard model (or its supersymmetric stringy extension), account for realistic large 
scale structure formations, as well as describe correct temperature properties of the cosmic 
microwave background. 

In the framework of strings, the theory of quantum gravity is well defined [1] (at least 
in certain cases) and can be considered at finite temperature as well [2-4]. Thus, in the 
stringy framework it is natural to attempt to describe the cosmological evolution of our uni- 
verse. In the very early times and when the temperature approaches the so called Hagedorn 
temperature [2-5] , we are faced to non-trivial stringy singularities indicating a high temper- 
ature non-trivial phase transition. In the literature, there are many speculative proposals 
concerning the nature of this transition [2,4,6-8]. 

Conceptually, it is quite exotic to understand the early stringy phase utilizing the stan- 
dard field theoretical geometrical notions [8-11]. Indeed, following for instance Refs [9,11], 
the notions of geometry and even topology break down in this very early-time era. Only 
by considering stringy approaches based on 2-dimensional conformal field theory techniques, 
there is a hope to analyze further this phase. Actually, some conceptual progress in this 
direction has been achieved recently [11], but we are not yet in a position to give quantita- 
tive descriptions of our stringy universe at these very early times. This stringy-conceptual 
obstruction however is not an obstacle to the study of the cosmological evolution of our 
universe for much later times than the Hagedorn transition era, namely for late times such 
that i > ijf, [12-14]. 

A way to bypass the Hagedorn transition ambiguities consists in assuming the emergence 
of three large space-like directions for t ^ tn, describing the three dimensional space of our 
universe, and possibly some internal space directions of an intermediate size, characterizing 
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the scale of the spontaneous breaking of supersymmetry, via geometrical fluxes, or even other 
(gauge or else) symmetry breaking scales [15]. Within this assumptions, the ambiguities of 
the "Hagedorn transition exit at t^;" can be parametrized, for t > Ie, in terms of initial 
boundary condition data at Ie- This is precisely the scope of the present work, were we 
focus on the intermediate cosmological era: < t < namely, after the "Hagedorn 
transition exit" and before the electroweak symmetry breaking phase transition at tw 

Modulo the initial boundary condition data at Ie, the stringy description of the in- 
termediate cosmological era, turns out to be under control, at least for the string back- 
grounds where supersymmetry is spontaneously broken by "geometrical fluxes" [16, 17]. 
These fluxes are implemented by a "stringy Scherk-Schwarz mechanism" involving n-internal 
compactified dimensions which are coupled non-trivially to some of the R- (super-) symmetric 
charges [18-20] [8,13,14]. Furthermore, the finite temperature effects are also implemented 
by a stringy Scherk-Schwarz compactification attached to the Euclidean time circle of radius 
(3 = 1/T [2-4,8,13,14]. Here, the R-symmetry charge is just the helicity, so that space-time 
bosons are periodic while fermions are anti-periodic along the Euclidean time circle. 

The breaking of supersymmetry and the quantum canonical ensemble at temperature T 
give rise to a non-trivial free energy at the one-loop string level. The latter is evaluated in 
a large regularized 3-space volume V3. This amounts to considering the free energy of a gas 
of approximately free strings, 

1 /7 (string) 

where Z = TYe~^^ is the partition function at finite temperature and Z[''_l^'^p^ is the one- 
loop connected graph i.e. genus one vacuum to vacuum string amplitude computed in the 
Euclidean background. In these equalities, we use the fact that the free energy of an infinite 
number of degrees of freedom can be found by a single string loop computation [21] [12-14]. 
The pressure derived from the free energy, P = —dF/dV^, has been computed for heterotic 
and type II model^ Supersymmetry is spontaneously broken by introducing geometrical 
fluxes d la Scherk-Schwarz involving n internal directions and by finite temperature effects. 
The pressure acts as a source occurring at 1-loop that backreacts on the original classical 
space-time metric and moduli fields background. In some cases, the perturbed equations of 

^Type I models realized as orientifolds of type IIB can also be considered. They can be dual to heterotic 
ones with an initially A/" = 4, 2 or even 1 supersymmetry [22] . 



2 



motion for the metric and moduli are shown to admit a cosmological evolution describing 
a radiation-like dominated era [12-14, 23] . The latter is characterized by a time evolution 
where the temperature T{t), the supersymmetry breaking scale M{t) and the inverse scale 
factor l/a{t) remain proportional [12-14,23]: 



Very different time evolutions also exist when temperature effects are neglected [14]. In this 
case, the supersymmetry breaking scale evolves like: 



In this work we would like to examine the late time evolution of the universe in the 
intermediate cosmological era, ts < t < ty,, i'n terms of initial boundary conditions (IBC) 
after the "Hagedorn transition exit at ig" • We show that the Radiation Dominated Solution 
(RDS) and the Moduli Dominated Solution (MDS) obtained for particular IBC are actually 
generic. For instance, for n = 1, there is a whole basin of IBC such that the resulting 
cosmological cvohitions converge - depending on the model - either to the RDS or to the 
MDS. The first attractor corresponds to an expanding evolution controlled by a Friedmann 
equation H'^ oc 1/a^, where H is the Hubble parameter. The second one describes a Big 
Crunch behavior dominated by the non-thermal radiative corrections. In addition, all models 
admit a second basin of IBC that yields to another Big Crunch behavior, dominated by the 
moduli kinetic energy i.e. where oc Within the scope of our present work, the 

Big Crunch evolutions can be trusted as long as Hagedorn like singularities are not reached. 
Our analysis is also limited to the cases where the initial supersymmetry is higher or equal 
to two, M >2. Although the M = 1 case is quite similar in some respects, it is much more 
involved in issues like the appearance of an effective cosmological term. It will be analyzed 
elsewhere. 

Restricting to cases with J\f > 2 initial supersymmetry, the attraction to the radiation 
era depends on the IBC: The convergence to the RDS is either exponential or via damping 
oscillations. During this process, periods of accelerated expansion can occur. However, they 
are characterized by an e-fold number of order 0.2, which is too low to account for realistic 
astrophysical observations. It would be very interesting to extend our analysis to models 



Radiation Dominated Solution : T{t) oc M{t) oc l/a{t) 




(1.2) 



Moduli Dominated Solution : M{t) (xl/a{tf^^/'' with T{t) <^ M{t) and (x 
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with A/" = 1 initial supersymmetry and see if the e-fold number is larger enough for an 
inflation era to occur. 

In Section 2, the case of a supersymmetry breaking using a single internal direction {i.e. 
n = 1) is studied in details. We classify the models and find local basins of attraction. 
This analysis is completed by a more general exploration of the set of IBC by numerical 
simulations. In Section 3, the attraction mechanism to the RDS is extended to models 
where n = 2 internal directions participate in the supersymmetry breaking. We determine 
in Section 4 the e-fold number of the periods of acceleration. Our results, their limitations 
and our perspectives are summarized in Section 5. 

2 Susy breaking involving n = 1 internal dimension 

In this section, we summarize first the results and notations introduced in Ref. [13] and 
then we analyze the attractor mechanism. The starting background can be the heterotic, 
type I or type II superstring compactified on a 6-dimensional space Ai that can be either 
S\R4) X or S^Ri) x x K3. The sizes of and x K3 are chosen to be of order the 
Planck lengthj^ The supersymmetry breaking is introduced by "geometrical fluxes" via the 
coupling of the momentum and winding numbers of the Scherk-Schwarz circle {R4) to an 
R-symmetry charge [3,13,14,18-20]. The radius of the circle i?4 controls the supersymmetry 
breaking scale M. In all applications, we will replace K3 by its orbifold limit, K3 ~ T^/Z2. 

In the heterotic and type I cases initially A/" = 4 or A/" = 2, the geometrical fluxes 
break spontaneously all the supersymmetries. A/" — > 0. In type II, when the S^lR^) coupling 
involves non-trivially left- and right- moving R-symmetry charges, all initial A/" = 8, A/" = 4, 
or A/" = 2 supersymmetries are spontaneously broken. When the S'^(i?4) coupling is left- 
right- "asymmetric" i. e. involves non-trivially a left- (or right-) moving R-symmetry charge 
only, half of the supersymmetries still survive: A/" A/'/2. 

^We relax this hypothesis in Ref. [15]. 
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2.1 General setup, gravity- moduli equations 



In all cases the computation of the free energy involves for the heterotic or type II superstring 
the background: 

S\Ro) xT^x S\R4) X with either = or M5 = x /Z2 . (2.1) 

The factor stands for the external space with regularized volume V3 sent to infinity once 
the thermodynamical limit is taken. ^^(-Ro) is the Euclidean time circle of period (3 = 2ttRq 
that defines the inverse temperature (in the string frame). In this direction, the Scherk- 
Schwarz circle amounts to coupling the momentum and winding numbers of the Euclidean 
time circle to the space-time fermion number [2-4,8,13,14]. The computation of the free 
energy density is well defined if the moduli Rq and -R4 responsible for the spontaneous 
supersymmetry breaking are larger than the Hagedorn radius Rh-, which is of order 1 in 
string units. At the Hagedorn radius, some winding modes become tachyonic so that the 
free energy develops a severe IR divergence [2-4,8, 13, 14]. If our goal is not to describe 
the phase transition that is occurring at the Hagedorn radius, we can suppose that in the 
intermediate cosmological era, Rq and R^ are large. In Refs [13,14], it is shown that up to 
exponentially small terms of order e~'^^° or e"'^^"', the pressure derived from the free energy 
density takes the following form [12-14,23], 

P = T*p{z) where = ' (2.2) 

with the temperature T and supersymmetry breaking scale M defined as, 

T= ^-^=, M= ?^^=. (2.3) 

In these expressions, the factors involving the dilaton in 4 dimensions, ReS* = e~^'^°, arise 



when we work in Einstein frame. In (2.2), p is a linear sum of functions of z with integer 
coefficients ut and ny [13,14], 

p{z)=nTfiz)+nvfiz), (2.4) 

where 

/(^) = 4^ 5Z [(2mo + l)2e2^ + (2m4)2]5/2 ' = ^''/(-^) ' (^.5) 

mo, 1714 \ / J 



and the sums are over the integers in Z. tit is the number of massless states of the model. If 
the Scherk-Schwarz circle S^{R4) acts non-trivially on the left- and right- moving R-charges, 
nv satisfies the inequality: 

-1< — <1, (2.6) 

and depends on the choice of R-symmetries used to spontaneously break all the supersymme- 
tries [13,14]. For an asymmetric 5'^(-R4) coupling i.e. trivial on the left- (or right-) moving 
R-symmetry charges, the supersymmetry breaking is partial and thus, there is no contribu- 
tion / to p. In that case, M is the mass of the Af /2 massive gravitini. In the following, 
we treat both kinds of models keeping in mind that ny = for asymmetric R-symmetry 
coupling to the momenta and windings of the 4*^ internal direction. 

From a low energy effective field theory point of view, the pressure P is minus the effective 
potential at finite temperature. It is a source for the quantum fields, namely the dilaton (pD 
and the modulus R^, coupled to gravity [12-14,23]. It also introduces the "thermal field" Rq 
to which there is no associated quantum fluctuation, as can be seen from the string model 
in light-cone gauge where oscillator modes in the direction are gauged away. The source P 
is a 1-loop correction that backreacts on the classical 4-dimensional Lorentzian background, 
via the effective field theory action. 



S = di'x\J—g 



1^) uv ( 1 1 ,1 duRid„R4 



(2.7) 



where we write explicitly the fields with non-trivial backgrounds only. Note that computing 
the 1-loop correction to the kinetic terms is not needed. Indeed, the wave function renor- 
malization of the moduli, (to recast the kinetic terms at 1-loop in a canonical form), would 
introduce an additional correction to the pressure at second order only. Thus, the full low 
energy dynamics at first order in string perturbation theory is described by the above action. 



(2.7). 



Before writing the equations of motion, it is useful to redefine the fields, 

|(0fl-lni?4), 0x := ^(20D + lni?4), (2. 



so that the action becomes: 

S = d^Xy 



(2.9) 



6 



where the pressure and the supersymmetry breaking scale M take the form, 



(2.10) 



. M e""^ . , /3 

= — , M = — with a=\-. 
T ' 27r V 2 

This shows that in the 2-dimensional moduh space {(PdjRa), the effective potential lifts 

the "no-scale modulus" 0, while keeping flat the orthogonal direction 0^. Assuming an 

homogeneous and isotropic cosmological metric with a flat 3-dimensional subspace, 



ds"^ = -Nitfdt^ + a{tf [dxl + dxl + dxfj 
the Friedmann equation in the N{t) = 1 gauge takes the form, 



3H' 



^12 , lj,2 



2" 2"^ 

where the energy density p is defined by the identity 

dP 



+ P: 



p + P = T 



dT 



(2.11) 



(2.12) 



(2.13) 



and P is understood as a function of the independent variables T and M. As shown in 
Ref . [14] , the expression of p follows from the variational principle applied to a gauge choice of 
the function N = 1/T. It is fundamental to stress here that the energy density defined by the 
variational principle is identical to the one derived from the second law of thermodynamics. 



From Eqs (|2^ and ( |2.13[ ), it follows that: 

p = r(z) where r = 3p — Pz 



(2.14) 



where 2;-indices stand for derivatives with respect to z = ln(M/T). Varying the action with 
respect to the scale factor a or the scalar fields and (f)±, we obtain: 



2H + 3H^ 
+ 3H(p = 
0x + 3iy0j 



1 > 



dP 



a{3P - p) 



72 



(2.15) 
(2.16) 
(2.17) 



where c_l is an arbitrary integration constant. Note that the last identity in Eq. (2.16) is a 



consequence of the fact that P is a dimension four object, {MSm + TdrjP = 4P, and from 



the definitions (2.10) and (2.13). The conservation of the total energy-momentum tensor 



follows from the Einstein equations (2.12) and (2.15): 



P + P 



0. 



(2.18) 



Using further the equations of the fields and we obtain from Eq. (2.18) that the 
{p, P}-system is not isolated but is coupled non-trivially to the "no-scale modulus field", 
[24,25]: 

p + 3H{p + P) + a^{3P-p) = 0. (2.19) 



2.2 Classification of the models, effective potential approach 



Our aim is to analyze the time behavior of the system described above. For this purpose. 



one can choose four independent non-linear equations, for instance: Eqs (2.12), (2.16), (2.17) 



and (2.19), any other choice being equivalent. Furthermore, we find convenient to express 



the dependence in time in terms of a dependence in the logarithm of scale factor a{t): 

dy dy 



A := In a{t) 



^ dt dX 



(2.20) 



for any function y. In particular, the definition z = acj) — ln(27rT) implies T/T = H{a(f) — z) 
that can be used to write 



p = HT^ 



+ (r^ - Ar)°z ) . 



(2.21) 



This equality is useful to express the conservation of energy-momentum (2.19) into a relation 
giving in terms of z and I, 

4r — r. 



ad) = Az(z)z — 1 where Aziz) 



3(r + p) ' 



Integrating, one obtains 



M 



(2.22) 



(2.23) 



where A{z) involves an integration constant. Using the above relations, the Friedmann Eq. 



(2.12) becomes 



m' = lH^ [{Az{z)z - If - 1) + 



9 e4(-^(-)--) , , 9 ci 
+ 



or, 



r{z) 



3- }C{z,z,(j)^) 



with /C = ^(^,(^)|-l)' + l0l 



(2.24) 



(2.25) 
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The equations for and (2.16) and (2.17), take the form, 
r{z) 



3 - IC{z,z,(f)^) 
r{z) 



AJzyz+A.Mz') + "-^^^^^ 



AJz)z + VJz) 



<Pi H ^ 



(2.26) 







. 3 - }C{z,z,(f)_^) 

For the derivation of the above equations, we used the relation y = Hy + H'^°y for y 



or 



as well as the Eqs (2/22), {2/25) and the linear sum of {2A2.) and (2AI)). As in Ref 



we have introduced the notion of an effective potential V{z), defined by its derivative (and 
using = 3/2): 

V,iz) = riz) - 4p(z) . (2.27) 

We are going to see in the following section that any solution of the system with constant 
z, i.e. z = Zc where Zc is an extremum of V{z), is giving rise to a particular RDS. The 
existence of such a critical point depends drastically on the shape of the potential V{z), 
which is determined by its asymptotic behaviors for z ±oo. It is obvious from Eq. ( 2.4[ ) 



and the definition (2.13) of p that V{z) depends only on two parameters, namely, nx and 



ny. The asymptotic behavior z ±oo can be easily derived: 
• {nv/nr) 0, 



V{z) ~ -e 

2— » + 00 



4z 



riy 
riT 



riTcl 



with 



r(5/2) 



TT 



5/2 A^\2m + lV> 

m ' ' 



(2.28) 



{nv/nr) ^ -1/15 = - {^'mj^) / {^mj2^)^ 



V{z) 



Z^ — OO 



4 



, with C4 



r(5/2) 2 



(2.29) 



riT 15/V3""V' " 7r5/2 3 ^ (2m + 1)^ ' 

According to the ratio nyjiiT., three distinct cases can actually arise, (see Fig. [T]), 

• Case (a) : -1 < {nv/nr) < -1/15 

V{z) increases monotonically from to +00. 

• Case (h) : -1/15 < {nv/nr) < 

V{z) admits a unique minimum Zc, with p{zc) positive. When {ny/nT) 0_ the 
critical value Zc goes to +00; when (nv/nr) — > (—1/15)+ the critical value Zc — ^ —00. 
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• Case (c) : < (nv/nr) < 1 

V{z) decreases monotonically from to — oo. 



v 

Case (a) 




Z 







Case (b) 




Case (c) 



Figure 1: The shape of the potential V{z) in Eq. (2.26) depends on the ratio nyjnT. The Cases (a), (h) 
and (c) correspond to —1 < ny jnT < —1/15, —1/15 < ny/nT < and < ny/nx < 1, respectively. 



Qualitatively, one can expect that the system represented by z can slide along its potential. 
It may run away in Case (a)., z(t) — oo, be stabilized in Case (b), z{t) Zc, and run 
away in Case (c), z{t) +oo. Moreover, these behaviors should imply new ones by time 
reversal. In particular, other run away behaviors are expected: z{t) —>■ +oo in Case (a), 
z{t) ±oo in Case (b), and z{t) — oo in Case (c). However, the study of the run away 
2; ^ — oo is out of the scope of the present work. This is due to the fact that this limit 
amounts to R4 ^ Rq > 1, so that the thermal system should be studied in 5 dimensions. In 
the following, we will only consider the other behaviors, namely the stabilization of z and 
the run away z — > +00. Since the latter amounts to -Rq ^ -R4 > 1? it is consistent with our 
analysis in four dimensions but with negligible thermal effects. 



2.3 Local attraction to the RDS 



Clearly, the first equation in (2.26) admits a constant z solution whenever V has an ex- 



tremum, i.e. in Case (b). It follows from Eq. (2.23) that M, T and 1/a are propotionnal 
during the evolution. 



M{t) = e'^T{t) 



a{t) ■ 



(2.30) 



The time dependence of the scale factor is dictated by the Friedmann Eq. (2.24) that 
simplifies to: 



3H' 



a4 a6 



where c 



9 



9 



(2.31) 
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Integrating, t can be expressed in terms of the scale factor, 



t{a) = to 



V3 



ao 



x^ 



Va > 0, 



(2.32) 



up to an overall sign. There are thus two monotonic solutions mapped to one another by 
time reversal. For the expanding one, the universe life starts in a "moduli kinetic energy 
dominated era" characterized by a total energy density ~ c^/a^, followed by a "radiation 
dominated era" with total energy density ~ c^/a^. For large enough times, the evolution 
always enters in the radiation era: It is asymptotic to the critical solution z 







i.e. (2.30) where 



i{t) = Vix 



1/4 



0. 



(2.33) 



To study the stability of this solution, we analyze its behavior under small perturbations. 
In terms of the variable A = Ina, we define, 



z{X) = z, + e{X) 



for small e and e±. At first order, the system (2.26) becomes 

P-Pzz 



oo 
6 



oo o 



^6 = 

0. 



where 



^ = 10 



I9p + pz 



(2.34) 



(2.35) 



The coefficient ^ is actually a function of the model dependent parameter ny/fiT, since the 
latter appears explicitly in the definition of p and implicitly via Zc- Numerically, one finds 
that ^ is always positive. It increases from to +oo when nv/nx varies from —1/15 to 0. 



The solution of the system (2.35) is 

r . .-A(i+v'i^)/2 



ci e 



+ C2 e 



-A(l-v'T=45)/2 



^±(A) = cse ^ , e{X) = < 



-A/2, 



-A/2 



Ci + C2 A) 



if < ^ < 1/4 
ife = l/4 



(cicos(Av/e^I74) +C2sin(Av/e^T74)) ife>l/4 

(2.36) 

where ci^2,3 are integration constants determined by the IBC. Since e{X) and e±{X) converge 
to for large A in all cases, the critical solution is stable under small perturbations. In 
other words, the solution arising for arbitrary IBC such that ^(0), e{0) and e±{0) are small 
is attracted toward the critical one, z = z^, <p± = 0. The behavior of e{X) is oscillating with 
damping when < ^ < 1/4, and exponentially convergent when ^ > 1/4. The special value 
^ = 1/4 corresponds to ny/riT — —0.0246. 
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2.4 Global attraction to the RDS 



Since we have shown the existence of a local basin of attraction around the critical solution 



given in equations (2.30) and (2.33), we want to know if this property is valid for more 

Q o 

generic IBC. Actually, any initial values (^o, zq, 0_i_o) at A = are allowed, as long as the 



positivity of the right hand side of the Friedmann Eq. (2.25) is guarantee^ i.e. 



rizo) 



> 0. 



3 - /C(zo,^o,0±o) 



(2.37) 



We have studied numerically the system (2.26) for generic IBC satisfying the previous bound. 
Implementing in the code a positive increment for A = Ina simulates the phases of the 



evolutions where the scale factor increases. We observe that the right hand side of (2.25) 



never changes of sign. This means that the scale factor a{t) is monotonic]^ We always find 

o 

that the solutions satisfy z{X) — » Zc and 0_l(A) ^ as A — > +oo, after possible oscillations. 
To illustrate the convergent behavior for large oscillations, we present a numerical example 
on Fig. |2} It is obtained for ny/riT = —0.02 that corresponds to Zc — 0.272, with IC 

Q o 

{zq, zo, (p±o) = (0.4,0.8,0.5). We conclude that for generic IBC, the expanding evolutions 




2 4 6 8 10 12 

A 

o 

Figure 2: Examples of damping oscillations of z{\) (solid curve) and convergence to zero of (dotted 
curve) illustrating the dynamical attraction toward the critical solution z = Zc — 0.272, 0_l = 0, obtained for 

nv/riT = —0.02. The initial conditions are {zq, Zo, 4>^q) = (0.4,0.8,0.5). 



are attracted by the RDS that satisfies z = Zc, (t)± = i.e. 3H^ = c^./a^ in Eq. (2.30) 



^IBC such that the right hand side of (2.25) is negative correspond to solutions in Euchdean time. 
^In general, a change of sign in the Friedmann equation during a simulation would be a numerical artifact 
meaning that H actually vanishes when the scale factor reaches an extremum. 
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The local behavior (2.36) also shows that the trajectory z = Zc, (px = is repulsive, when 
A decreases i. e. for evolutions where the scale factor is contracting. The Big Crunch solution 



obtained by time reversal on (2.33) is thus only formal, since it is highly unstable under small 
fluctuations. What we observe by numerical simulations is that for generic IBC, when we 
choose a negative increment for A, then either z +oo, or z oscillates with a greater and 
greater amplitude. To understand better the dynamics of the contracting evolutions, we are 
thus led to analyze carefully the system in the regime z ^ 1. 



2.5 Attraction to Big Crunch eras 



We have seen that when the potential (2.27) admits a minimum z^. {Case (b)), the expanding 
cosmological evolutions for arbitrary IBC satisfy z{t) Zc- The shape of the potential in 
Case (c) of Fig. [l] suggests that when ny > 0, the dynamics could admit a run away 
behavior z{t) +oo. We also noticed at the end of the previous subsection that the 
contracting evolutions in Case (h) involve the large and positive value regime of z. For these 
reasons, we need to study the system for z ^ 1. In this limit the thermal effects 0{T^) 
are sub-dominant compare to the supersymmetry breaking effects (9(M^). The expansion 
oi p{z) contains two monomials in e^, 

r(5/2) 2>^' 1 



p{z) = nye ""cg + {utcI + nycl) + ■■■ where cl 



5/2 3^ (2m) 



4 ' 



(2.38) 



and defined in Eqs (^2.28| and (^2.29^. The dots stand for exponentially suppressed 

terms 0(e~^"^). It follows that Az{z) ~ 1 and the conservation of the energy- momentum Eq. 



(2.22) becomes 



T 

z = M/M + l ^ f^~^ ciT = aoTo, (2.39) 

where QoTq is a positive constant depending on the IBC. From the expansion of p{z) for 



large z in Eq. (2.38), it is clear that the behavior of the system is drastically different when 
ny 7^ and when ny = 0. 

• ny 7^ 0, z ^ 1 : 



We have p ~ — P ~ — M nyCg. The Friedmann Eq. (2.12) and the linear sum of Eqs (2.12) 
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and (2.15), together with the matter field Eqs (2.16) and (2.17), become 



a(f) + 3Ha<p = Aa^M^nycl , 



+ ?,H(t)i_ = 



(2.40) 
(2.41) 
(2.42) 



The above equations would also arise if we had considered the system at zero temperature. 
This is due to the fact that when temperature effects are not taken into account, the sources of 
Einstein gravity always satisfy p = —P equal to minus the effective potential (see Eq. (2.13 )). 
However, this does not mean that in our case we shall find that T(t) 0. Actually, the 
regime z ^ 1 i.e. M ^ T corresponds to thermal effects screened by radiative corrections. 



even if the temperature can be large. Combining Eq. (2.41) and Eq. (2.42) one determines 



the scalar fields in terms of the scale factor, (with a = -\/3/2), 



-4a'H + 



M 



ft dt' 



a-" 



a 



6 



and 



V2 



a3 



(2.43) 



M^nvc° . 



(2.44) 



where c^p and c± are constants depending on the IBC. The Friedmann Eq. (2.40) takes then 
the form 

3 V a^J 

We will consider separately the cases = and 7^ 0. 
(i) C0 = : 

M = Mo 



3i7^ = -4 + % with c.-^ 



My,al' , (2.45) 



aJ ' a" ' a'"' 3 

where M^a^ is a positive constant. This solution exists only for ny > . If we also have 
c_L = 0, up to time reversal, the evolution describes a Big Crunch occurring at a finite time 

^BC5 

a{t) = ao (tec - ty^'^ X (16Mo^nyc°)i/2^ , t < tec • (2.46) 

Since a ~ 1/T, it follows that oc aM oc 1/a^, which is large if t ^ tec and implies a run 
away of z, consistently with our hypothesis z ^ 1. However, close to the Big Crunch, the 
scale M (as well as T, even if T -C M) is formally diverging. This implies that the present 
analysis has to be restricted to the regime where M < Mfj i.e. z < {5/Q)ln M^, where 
Mh is the supersymmetry breaking scale above which a Hagedorn-like transition occurs. To 
bypass this limit, one could try to extend our study to Hagedorn singularity free models [8]. 
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For c_L 7^ 0, the time can be expressed as a function of the scale factor, 



t{a) 



A8 



where U 



\cm\ 



,12 



Va < a„ 



Cm 



1/9 



(2.47) 



Formally, this solution describes a Big Bang at t = — that initiates an expansion era. 
The latter stops when the scale factor reaches its maximum a,^^^^ at t = 0. The universe then 
contracts till t = tj^c, where a Big Crunch occurs. As before, the last era of this evolution 



can be trusted for t ^ t^c, where the solution (2.47) is asymptotic to (2.46). Thus, we have 
shown that for IBC such that zq ^ 1 and = 0, the dynamics is attracted to the MDS 
described by the Big Crunch solution of Eq. (2.46) with = c_l = and 3H^ = cm/o?'^- 

(ii) C0 7^ : 

We consider IBC such that zq 3> 1, with generic ^ 0. The sign of nv 7^ is arbitrary and 



we want again to find the basins of attraction of the dynamics. Instead of Eq. (2.40), it is 



more convenient to use the combination of Eq. (2.40) and Eq. (2.41), 

1 > 



2^ 



2 



that becomes, utilizing the expressions (2.43) for and 

1 



E 



Introducing the quantity a. 



4q 



Eq. (2.49) takes the form: 



6H 



3 

4ca 



4:Ca 



da 



-9 



da 



a2-a+^(l + 3(cx/c<^)2) a 
valid in the z ^ 1 limit. In the same limit the field equations become: 

^ -(M/M)°-M/M -6 = 



(2.48) 



(2.49) 



(2.50) 



(2.51) 



|(M/M)2 + l02_3 



(2.52) 



where the signs of ny and the fractions in front of the second derivatives must be the same. 



so that the right hand side of Eq. (2.25) is positive. This condition is required to study the 
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evolutions in real time (and not Euclidean). This differential system is easily shown to be 
equivalent to 

' ?,{M/My = {M/M + 6)V{M/M) 

° , , (2.53) 
0^ = — (M/M + 6) ^ ^ 

where c±/c^ is an arbitrary constant denoted this way to make contact with the conventions 



of Eq. (2.43), and the polynomial 



V{x) = [1 + 3{c±/c^y] + 36 {c±/c^y x + 9 [12{c±/c^y - 1] is of the sign of nv- (2.54) 



Using the first equation in (2.43), one obtains 



M/M + 6 



(2.55) 



that can be used to identify consistently (2.51) with the first equation of (2.52). By the way. 



the real time condition translates to the condition that ny and the denominator of the left 



hand side of (2.51) have a common sign. We are going to solve this equation separately for 



ny > and ny < 0. 
- For ny > 



i) When (c_l/C(^)^ > 1/9, the denominator of (2.51) has no real root. Integrating, one 
obtains 



a 
ao 



exp 



9^9(ci/c^)2-l 



tan^ 



-1 



40--2 



a 



16* 



1/18 



(2.56) 



that is used to draw a versus a on Fig. |3]a. Suppose > 0. If at some time a > 0, 



the definition of a in Eq. (2.50) implies that a increases. Thus, the scale factor reaches 
a maximum a^^^ where a vanishes. We then pass into the a < lower half plane where 
a decreases and converges to 0. Since a —>■ — oo and a 0, one has a —>■ — oo, which is 
interpreted as a Big Crunch occurring at a finite time tec- Similar arguments apply when 
< : If cr < 0, the scale factor increases up to a^^^ before converging to 0. This 
corresponds to a Big Crunch occurring at some finite time tec- To summarize, when the 
IBC are such that > (c^ < 0), the system follows the curve on Fig. [3]a from up to down 
(down to up). The universe starts from a Big Bang, reaches a maximum size, and ends with 
a Big Crunch. Of course, only the parts of this evolution consistent with the hypothesis 



16 



a) G 




max 



a 




max 



a 



Figure 3: a, defined in Eq. (2.50) versus the scale factor a. Fig. a) corresponds to {c_i_lc^)'^ > 1/9 md 
Fig. h) corresponds to [c^jc^Y — 1/9- The solid lines are associated to the case ny > 0, while the dashed 
one is associated to the case uy < 0. 



2; S> 1 can be trusted. In particular, when we approach tsc, one has a 
s = sign(c<^), and Eq. (2.56) imphes, 



a(t) ^ (tec - t f'' X ((16/3) |c^| al 6^^^/2)1/12 



-soo, where 



(2.57) 



This result can be used to show that the integral in M in Eqs (2.43) is bounded when t t^Q 
so that a consistent run away oc 1/a^ — > +oo is found. An important remark then follows. 



The behavior (2.57) implies that c^/a^ is dominated by H in Eq. (2.43), showing that the 
dynamics for {c^/c^Y > is always attracted by the solution for c,^ = c_l = 0. 



ii) When (c_l/c0)2 < 1/9, the left hand side of Eq. (2.51) has two real roots. Integrating, 
one has 



for 



for 







a 












ao 






a 






ao 



0"_ 



1 



I 9{<T+-<T_) 



g 18(<T-l/2) 

Iff - l|l/9 
1"^ 2 1 



1 



2 4 



l-9(c^/c^)2, (2.58) 
(2.59) 



where cr_ < a < a^. Fig. [3^ shows cr as a function of When c,^ > 0, the arguments used 
in case i) apply on the branch cr < o"_, where the system evolves from up to down. The 
scale factor grows up to a^-^^^ before converging to 0. Similarly, when < and the system 
follows the branch a > 0"+, the scale factor converges to 0. In both Big Crunch 



^To be precise, the slope at the point {a, a) — (0,cr_) is infinite when (1 — 4/361)/9 < (cj_/c0)^ < 1/9, 
finite when {c^/c^f = (1 - 4/361)/9 and vanishing when {c^/c^f < (1 - 4/361)/9. 
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occurs at tac and the behavior (2.57), with s = 0, is vahd. The dynamics is again attracted 
by the solution associated to = c± = 0. 

Some new features occur for < 0, when a < a-, (see Fig. The scale factor 

increases up to its maximum a,,^^^ before converging to 0, while a — cr_, implying a — > — oo. 
A Big Crunch at a finite time tsc occurs but since a cr_, one has 

a(t) ~ (tec-t)'/' X ((4/3)|c<^|(T_)^/^ t<Uc. (2.60) 

_3 5 

However, this behavior can be trusted as long as oc aM oc Itec ~ tl*"^- ^ — > +oo. This 
condition implies that the scaling (tec ~ ^Y^^ exists when 

1 /c \ ^ 1 

c^<0 and - (1 - 1/25) < j < - . (2.61) 

In this case, the evolution is not attracted by the solution c<^ = c_l = 0. Actually, Eq. (2.43) 
can be rewritten as 



a0 ~ — 2 In Itac — t\ + c—> +00 when t ^ t^^ ^ (2.62) 

so that the integration constant c is negligible. There is thus a new basin of attraction, to 



the solution (2.60) with c = i.e. such that 

^H^-"^ where c^ = ^c5a!, c = 0. (2.63) 

As for the previous Big Crunch behavior, the divergence of M (and T) is only formal since 
the present analysis supposes M < M^, where a Hagedorn-like transition occurs. On the 



contrary, for (c_l/c^)^ smaller than the lower bound of the range (2.61), we formally have 

3 5 

oc |tBc —tl"^"^- ^ — s> 0, which shows that the system actually exits from the regime z ^ 1 

we started with. Similarly, if > and the system evolves along the branch a > of Fig. 

one has a +oo and a cr+. This implies a ^ t^^'^ x ((4/3)C(^cr+)-'^/'^ and t +oo. 

_3 5 

However, for large positive z, this implies oc aM (x. 1'^'"+ ^ — > 0, which shows that the 
system actually exits in the regime z ^ 1. 

Summary for ny > ^ '■ 

We have shown analytically that when 1, either z —>■ +oo or z quits the regime z ^ 1. 



To study the dynamics at finite z, we use again a code that implements Eqs (2.26). The 
phases of expansion are simulated when a positive increment for A is chosen. For generic 
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IBC, we observe that the right hand side of the Friedmann equation (2.25) changes sign at 
some finite a = cimax- This actually means that the scale factor has reached a maximum, (the 
Hubble parameter H vanishes), and that the evolution enters into a phase of contraction. 
Such eras are simulated by choosing a negative increment for A and we observe that for 
generic IBC, z ends by running away. As a conclusion, the generic IBC admit two basins of 
attraction to MDS. The evolutions converge either to the solution where = c_l = i.e. 
such that 3H^ = c\^/a^^, or to the solution with c = and 3H^ = c^/oP- 

- For < ; 

For IBC such that ^ the real time condition translates to cr_ < cr < cr+, (see Fig. ^p). 
If > 0, the definition of a in Eq. (2.51) reaches a > 0. It follows that a +oo and 

^3 5 

cr 0-+. Thus, a ^ t^/^ X {{A/3)c^a+y'^ where t +cx), and oc aM oc ^ 0. 
We conclude that the system exits the regime z ^ 1. When —1/15 < ny/nx < {Case 



(b)), this result is compatible with Sect. 2.3, where we found that the expanding evolutions 
satisfy z ^ Zc- 

If on the contrary < 0, one has a < 0. It follows from Fig. [3^ that a —>■ and a ^ a_, 
so that d — oo. The dynamics is attracted by the Big Crunch solution (2.60) i.e. ( |2.63[ ) 
that can be trusted when the conditions (2.61) are satisfied. For smaller values of (c^/c^)^, 
the system exits the regime z ^ 1. These conclusions are compatible with the numerical 



simulations of the contracting evolutions described at the end of Sect. 2A, for generic IBC 
in Case (b). We found there that either z — > +oo, as expected from the existence of the 



Big Crunch solution (2.60), or z oscillates with a larger and larger amplitude, implying that 
when z is large and positive, it can exit the regime z ^ 1. 

• ny = 0, z ^ 1 : 



We argued in Section 2.3 that when —1/15 < ny/riT < 0, the expanding evolutions are such 
that z is attracted by a critical value Zf.. Since Zf. +oo when ny/riT —>■ 0_, we expect 
the expanding solutions for = to satisfy z{t) +oo. Actually, this run away may be 
natural since the potential in Eq. (2.27) decreases linearly for large positive z, V ^ —zutcI. 



In fact, a numerical study of the differential system (2.26) with ny = confirms that for 
generic IBC, z enters the regime z ^ 1 when we simulate the expanding solutions. On the 
other hand, for the evolutions where the scale factor decreases, we find either z +oo or 
z oscillates with a larger and larger amplitude. In all these cases, we are thus invited to 
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consider analytically the system in the regime z ^ 1. 



In this limit, the thermal effective potential for in Eq. (2.16) vanishes (up to exponen- 
tially small terms in z). Both and (f)± are flat directions, and their kinetic energies involve 
two constants c^, c±, determined by the IBC, 



a I 



C4, 



V2 



(2.64) 



The energy density in the Friedmann Eq. (2.12) is p ^ 3T nxc^ and, using (2.39), one has 



3H' 



+ 



where = 371^04(00^0 



(2.65) 



This equation admits two monotonic solutions (see Eq. (2.32)), mapped to one another 
under time reversal. However, one can only trust them as long as z ^ 1. 



- For the expanding evolutions, since a(t) — > +00 when t +00, Eq. (2.64) implies 
that M{t) converges to a constant Mq. It follows that = M/T oc a +00, showing 
the validity of the solution for large enough time. We conclude that for generic IBC, the 



expanding evolutions are attracted to the RDS, Eq. (2.33), with M = Mq. 
- For the contracting evolutions, depending on and c±, we have 

a{t) ^ (tec - ty^^ X {3ciy/' , M ~ Mo 



t 

^BC 



"[l+3{c^/c<^)2]l/2 



t<tBc, (2.66) 



if the constraints 



are satisfied, (so that z 



C0 > 
+00 when t 



and 



C_L 



< 



(2.67) 



/4 \ ' 

j(t) ~ Vt;^x ( -cM , M~Mo 



Alternatively, if = c± = 0, the solution 

1/4 



(2.68) 



is also consistent with a run away of z. For all other values of c<^ and c±, z exits the 
regime z ^ 1. Our conclusions for the contracting evolutions are that there are two sets 
of generic IBC. The first gives rise to a run away of z that corresponds to the Big Crunch 



solutions (2.66) or (2.68). The second yields an oscillation regime in z, with larger and larger 
amplitude. 
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3 Susy breaking involving n = 



2 internal directions 



In this section we extend our analysis to models where the supersymmetry breaking is in- 
duced by = 2 Scherk-Schwarz circles 5'^(-R4) x S^i^R^) in the 4*^ and internal dimensions 
as in Ref. [14]. Thus, we are considering the heterotic or type II superstrings on 

S\Ro) X X S\R4) X S\R5) x , (3.1) 

where Ai^ is either or K3 ~ T^/Z2. As shown in [14], such models allow cosmological 
evolutions that describe radiation dominated eras, with stabilized complex structure modulus 
R5/R4 [14]. The present analysis could certainly be extended to the other classes of models 
described in [14], where asymmetric Scherk-Schwarz compactifications in the directions 4 
and/or 5 are considered. These models would involve run away solutions in the spirit of the 
Cases (c) and (a) for n = 1. 

Rq , R4 , -R5 are restricted to be large enough in order to avoid Hagedorn-like phase 
transitions. In this regime, the free energy and pressure are computed in Ref. [14] : 

P = T^p(z, Z) where = — , T = ^-^= , M = ^ , , = — . 

T 27tRoVR^ 27ry^Rjr5VR^ Ri 

(3.2) 

T is the temperature in the Einstein frame, M is the supersymmetry breaking scale which 
is given in terms of the geometric mean of the moduli i?4 and -R5, and Z parametrizes the 
complex structure modulus R5/R4. The quantity p{z, Z) is a linear sum of functions, with 
model-dependent integer coefficients, 

p{z, Z) = riTPiooiz, Z) + nQiopoiQ{z,Z) + nooiPooi(^, Z) + niiipiii(z, Z) , (3.3) 

where ut is again the number of massless states and 
2 v-^ e^" 

= ^ [(2mo + ~9o?e^^ + {2m, + ~g,fe-^ + {2m, + -g.fe^f ' ^^'^^ 

3.1 Gravitational and moduli equations 



The low energy effective action (2.7) with the kinetic term of R, taken into account can be 



written in terms of redefined fields as : 



S = d xv— ^ 



lR-l({d<pr+{d<p^r+l{dzr]+p 
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(3.5) 



where 



+ In Vi?4i?5 , P = M^e-^'p{z,Z) , M=^. (3.6) 



For evolutions satisfying the metric ansatz (2.11), the Friedmann equation is 
3H'^ = + 1^2 ^ 1 + r where r(2, Z) = 3p - p, , 
while the equation for the scale factor, 

2 4 2 



(3.7) 



(3. 



can be replaced by the conservation of the energy-momentum in a form similar to Eq. (2.19 ) 



p + 3H{p + P) + ^T\3p - r) + ZT^ pz = Q . 



(3.9) 



For the scalar fields, one has 



+ = (3p - r) 
Z + 3HZ = 2pz , 
0x + 3H^^ = 



(3.10) 
(3.11) 
(3.12) 



We proceed by introducing the variable A of (2.20) as in Section 2.2 to express 



A{z,Z)l + B{z,Z)Z -1 where A{z,Z) = ^ — ^, B{z,Z) = - ^f (3.13) 



3(r +p) 



3(r +p) 



that reaches, by integration, 

M=- 



a 



where := ^(z, Z)^ + B{z, Z)Z . 



(3.14) 



This result can be used to rewrite the Friedmann Eq. (3.7) in either of the following forms, 

2 g ^i[H^na) 



3H^ = Ih^ ({Xz + BZ-lf + \z-l 



a'* 



r{z,Z) + 



6 cj 
5 a 



6 



(3.15) 



or 



r{z,Z) 



o ° ° 

3 — }C{z, z, Z, Z, 0j_) 



where /C = ^ (^A{z, Z)z + B{z, Z)Z - Ij + ^ + i ^ 

(3.16) 
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while the equations for the scalars become 



o o 



A°z+BZ + A°z^ + {Az + BSzZ + BzZ^ \ + + V^^'^ = 



< — ^ z + ^ Z + V^/^ = 

3-/C 2 ^ 



' + — 



I 3-/C 



where we have defined 



Vj^\z,Z) = '^^ and \/f)(z,Z) = -2p^ 



(3.17) 



(3.18) 



3.2 Attraction to the RDS 



A wide range of behaviors can emerge from this system of differential equations. However, in 



the spirit of Section 2^, we choose to look for cosmological solutions with stabilized complex 
structures = Ro/V RiR^ and = R5/R4 i.e. satisfying {z,Z) = {zc,Zc). From the 



system (3.17), such solutions exist if the following conditions are simultaneously satisfied. 



v}'\z,,Z,) 

\^c, Zc, 



I.e. 2p{zc, Zc) +Pz{zc, Zc) = 
i.e. pz{zc,Zc) = 0. 



(3.19) 
(3.20) 



These equations define two sets of curves in the {z, Z)-plane, whose intersections determine 
the critical points {zc, Zc). Supposing such a solution exists, the corresponding cosmological 



evolution is characterized by JF = in Eq. (3.14) and thus 



M{t) = e'^ T{t) 



a{t) 



where 



= est. 



together with the Friedmann Eq. (3.15) 

^2 



3H' 



'^r _|_ 

a4 ^ a6 



R^it) = e^^ R,{t) 



6 



(3.21) 



where 4 = 6 e^(^-^^) p{zc, Zc) , d = - 4 



(3.22) 



The latter is qualitatively identical to the one found for n = 1 in Case (b). The evolution 
with expanding scale factor is asymptotic to the critical solution {z, Z) = {zc, Zc), Cm = 0, 



where a{t) is given in (2.33). 
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To study the local stability around a critical solution with Cm = 0, we define small 
perturbations around it, 



^(A) = z, + £(A) , Z(A) = Z, + E(A) , 0^(A) = 



(3.23) 



and consider the system (3.17) at first order in e, E and e±. The latter can be brought into 
the form 



oo o 



9p{'ip-pzz)+2p 



e + e + ^e + uoE = Q where ^ 



i + E + SE + l]£ = where S = - 2^ 

p 



2p{26p+pzz) 

iZc,Zc) 



PzZ {2pzz-'3p) 
2p (26p+pzz) 



P 



uu u 

, e±+e± = 0. 



(3.24) 

To proceed, we specialize on the set of models whose R-symmetry charges coupled to the 
winding and momentum numbers in the directions 4 and 5 are identical. For this set, it is 



shown in [14] that the function p in Eq. (3.3) takes the form 



p{z, Z) = utIpioo + Pni] + '^vboio + Pooi] , -1 < — < 1 



(3.25) 



and that the range of ny/riT can be divided in four phases whose frontiers are given by 
{-0.215, -1/31, 0} [14]. 



For ny/riT ^ —0.215, the constraint (3.19) has no solution 



For —0.215 ^ ny/riT < —1/31, the constraint (3.19) defines a non trivial curve, while 



the condition (3.20) is satisfied on the axis Z = 0. The two loci intersect at a point 
(^c,0). 



For —1/31 < ny/riT < 0, the constraint (3.19) defines a non trivial curve, while the 



condition (3.20) corresponds to 3 curves, one of which being the axis Z = 0. There is 



a unique intersection point (-2c,0). 



For < ny/riTi the constraint (3.19) has no solution. 



Thus, a solution to the equations of motion with constant complex structures (zc, Zc) exists 
when —0.215 ^ ny/riT < 0. 
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The class of models we consider have a symmetry R4 
This implies Pzz{zc,0) = and 



i?5 that translates to Z ^ —Z. 



9 4p-p2 



2 26p + p,, 



0, H 



Vzz 



1] = 0. 



(3.26) 



As in the n = \ case, one finds numerically that ^ and S are always positive. They in- 
crease from to +00 when nx/nv varies from ~ —0.215 to 0, which shows that e(A), E{X) 
(and ^_l(A)) are always converging to zero. This means that in the class of models under 
investigation, when a solution (z^, Z^ exists, it is stable under small perturbations. 

I I o 00 

A numerical study of the exact system (3.17) shows that for generic IBC {zq, Zq, Zq, Zq, 0_|_o) 



such that the right hand side of Eq. (3.16) is positive, the expanding cosmological solutions 
satisfy {z{X), Z{X)) {zc, Zc) when A +00. We conclude that the dynamics of the ex- 
panding evolutions is attracted by the RDS described by the critical solution {z, Z) = {z^, Z^), 
^± = I.e. 3H^ = clla^ (Eq. E^). 



4 An inflation era ? 

We have described how cosmological evolutions can be determined by the R-charges of the 
massless spectrum. In particular, we found that the dynamics can be attracted to a radiation 
era, with stabilized complex structures. In standard cosmology, such an era follows a period 
of infiation introduced to explain the observed fiatness and homogeneity of our universe, 
or solve the topological defect problem. It is thus important to analyze if our attraction 
mechanism admits periods of accelerated cosmology. 

Using the Einstein equations, the infiation condition a > or if + H'^ > can be 
translated into an inequality between the thermal sources p, P and the scalar kinetic energies 
(take Z = for n = 1), 



Using the variable A and Eq. (2.25) (or (3.16)), we obtain 



1 ^ / 1 \ 1 



p--p.<-2{p--pA ^-^ < , (4.2) 

where there are two conditions in one: First, p — Pz/Q has to be negative, and second it has 
to be negative enough compared to the right hand side that depends on the kinetic energies. 
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For the classes of models we considered for n = 1 (see Eq. (2.4)) and n = 2 (see Eq. 



(3.25)), one can show that a necessary condition to have p — Pz/Q < is < 0. We thus 



carry on our discussion in this case. It then follows that the quantity p — Pz/S in (4.2) is 
always positive. Altogether, a necessary and sufficient condition for an accelerated cosmology 
is: ^ 

nv <0 and /C < J where I = . (4.3) 



p 



\Pz 



Actually, when deriving (4.3) from (4.2), one also finds that X < 3. However, this is an 



empty constraint since the function I satisfies T < 1. 
For a susy breaking with n = 1 



From (4.3), in order to have solutions, we need T(nv/nT, z) > 0. This condition defines a 



domain in the {ny /ut, z)-p\aii where, by definition, X = on the boundary. However, it 
happens that X ^ 1 in the interior, as can be seen on Fig. |4]a. We can thus estimate the 



a) 



b) 




Figure 4: For one modulus z (Fig. a)): In the domain 2 > of the (ny /tit, z) -plane < ny/nT < 

0), X can he approximated by 1, its value for z 3> 1. For two moduli {z,Z) (Fig. b)): For any fixed 
—0.215 ^ ny jriT < 0, in the domain X > of the {z, Z)-plane, X can be approximated by 1, its asymptotic 
value in the sectors II, III , IV, V defined in Eq. (^.9). Fig. b), is drawn for ny /nx — —1/63. 



e-fold number by considering z large and positive. In particular, when —1/15 < riy/'^T < 0, 
the z's allowing periods of accelerated cosmology are always larger than Zc- When z ^ 1, 



the condition (4.3) simplifies to 

1 



{M/Mf + {c^/c^)\M/M + 6)2 < 1 



(4.4) 



where we have used the second equation of (2.53). An era of acceleration exists if {c^/c^f < 
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1/33. It begins and ends at most when MjM saturates the inequahty (4.4) i.e. equals 

-18(cx/c<^)2 ± v/3[l - 33(cx/c^)2 



6,e 



1 + 3(cx/c^)^ 



(4.5) 



To find the scale factors and Oe when this arises, we integrate the first equation of (2.53), 
= A- Ao = 7ln|M/M + 6| +7+ln|M/M-/+| +7^ ln|M/M-/_|, 



In 

where 



(4.6) 



7 



l + 3(cx/c<^)2 (/+ + 6)(/_ +6) 



-18(cx/c0)^ ± 3v/l - 9(cx/c^)- 

1 + 3(cx/c<^)2 
L 3 

7± = 



(4.7) 



l + 3(cx/c^)2 (Zi + 6)(/+-/_) 



Inserting the particular values (|4.5|) in (|4.6|), the maximum e-fold number is found to be 

If, — /_|_ 



e(c_L/c0) = In 



Ob 



7 In 



/e + 6 



4 + 6 



+ 7+ In 



Ih — 



7_ In 



le-L 



h-L 



that satisfies e{c^/c^) < e(0) ~ 0.227. Thus, even if periods of accelerated scale factor are 
allowed, they do not give rise to infiationary eras. 

For a susy breaking with n = 2 

For any fixed ny /ut < 0, the condition I{nv/nT, z,Z) > defines a domain in the {z, Z)- 
plan, as shown on Fig. |4^. Using the notations of [14], this zone spans four asymptotic 
sectors defined as follows, 

// : Z ^ —oo, z ~ rjZ, \ri\ < ^ , 
IV : Z ^rjz ^ +00, < ?7 < 2 , 



/// : Z rjz ^ -oo, -2 < r/ < 
V : Z ^ +00, z ~ -qZ, |?7| < ^ , 



(4.9) 



inside of which X ^ 1. We will evaluate the e-fold number by approximating p{z, Z) by its 
asjTiiptotic expansions, which are well defined in each sectors [14]: 



pii ^ e^'e-^^nvSl + e"'e~^'^{nTSl + nySl) + e^'e^^'^^riT + nv){Sl + ^4^) + 



III 

IV 



nvS^ + e^'e^^nviS^ + SI) + 2{nTS2 + ny^|) + 
e^^e^^ny^" + e'^'e-^^nyiSI + SI) + 2{nTS2 + nySl) + 



e e 



P 
P 

pV ^ e^'e^^nvS^ + e-'e^/^{nTSI + nySI) + e^'e'^^/^{nT + nv){S2 + SI) + ■ ■ 



-fO,e 
'4,5,( 



where the dots stand for exponentially suppressed terms and the coefficients SI 
stant£l 

^Their precise values are not needed in the following but can be found in the Appendix of [14] 



(4.10) 



are con- 
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It is interesting to note that the behavior of the temperature depends drastically on the 
location of the representative point of the system in the {z, Z)-plan. In sectors /// and 
IV, one finds ^ ~ 1, i3 ~ as in the 1 modulus case, while in sectors // (or V^), one has 
^~4/3, 1/6 (or —1/6). Using the energy-momentum conservation in the form of Eq. 



Ill U IV : aT = est. 



V : — — - = est. . 



(4.11) 



(3.13), this reaches 

: ^g^g -est., . .... , . . ^g^g 

However, the Eqs of motion for the scalar fields are identical in all these asymptotic sectors, 

1 



-(M/M)°-M/M -4 = 



i(M/M)2 + iZ2 + l02_3 

1 OO o 

— , Z -Z + 6 = 



|(M/M)2 + 1 Z2 + 1 02 _ 3 



(4.12) 







t l(M/M)2 + iZ2 + l02_3 
where the fractions in front of the second derivatives must have the sign of ny i. e. negative. 

o 

When M/M + 4 7^ 0, it is easy to show that there exist constants 7^ 0, c^, c± such that 



2(M/M)° = (M/M + A)V{M/M) 
Z-6 = 2 — (M/M + 4) 

0^ = ^2 — (M/M + 4) 
C4, 



(4.13) 



where 



V{x) = [1 + 2{c,/c^f + 2{c^/c^f] x2 + A[A{c,/c^f + 3cjc^ + A{c^/c^f] x 
+4[8(cz/c0)2 + I2cz/c^ + 8(c_l/c0)^ + 3] is of the sign of ny. 



(4.14) 



The condition (4.3) is also identical in all asymptotic sectors II,...,V, 

/C ~ ^ (M/M)2 + [{cz/c^){M/M) + 3]2 + {c^/c^y{M/M + Af < 1 



(4.15) 



where we have used the second and third equations of (4.13). A period of accelerated 
cosmology exists if 



14 



< Cz/ctf, < c+ where c± 



6±2v/2[l - 14(cx/c<^)2] 



(4.16) 
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It begins and ends at most when M /M saturates the inequahty (4.15) i.e. equals 



To find the corresponding scale factors Of,, Og, we integrate the first equation of (|4.13|) 
In 



where 



A- Ao = 7ln|M/M + 4| +7+ln|M/M-/+| +7_ ln|M/M-L 



-2[A{c,/c^f + 3c Jc^ + 4(cx/c^)'] ± ^1^{C+ ~ c,/ c^){c,/ ~ C.) 



with C± 



1 + 2{c^/c^f + 2(cx/c<^)2 
-6 ± v/7[3 - 10(cx/c<^)2] 



7 



7± 



1 + 2(cjc^)2 + 2(cx/c^)2 (Z+ + 4)(/_ + 4) 
2 1 



■ 1 + 2(c,/c^)2 + 2(cx/c<^)2 (/± + 4)(/+ - /_) 



Inserting (4.17) in (4.18), the maximum e-fold number is found to be 

^e + 4 



e{cz/ C0, c_l/ c^) = In — = 7 In 



^b + 4 



+ 7+ In 



4 — ^+ 



+ 7_ In 



L - L 



(4.18) 



(4.19) 



(4.20) 



that satisfies e{cz/c^,c±/c^) < e(— 3/2,0) ~ 0.189. The order of magnitude of the e-fold 
number for n = 1 and = 2 is thus the same. They cannot account for the astrophysical 
observations. 



5 Summary, limitations and perspectives 

The present work focuses on an intermediate cosmological era, Ie < t < t.^. The time tE 
corresponds to the exit of an Hagedorn phase i.e. after the notion of topology and geometry 
emerge, with 3 large space-like dimensions and possibly internal directions of intermediate 
scales. The characteristic time is associated to the electroweak phase transition, occurring 
before large scale structure formations in the universe. 

To describe the intermediate era, we consider fiat classical 4-dimensional space-times 
within the context of superstring theory, where A/" = 8, 4, 2 supersymmetry is spontaneously 
broken. Our aim is to determine the different cosmological behaviors that can be induced 
when three kinds of effects are taken into account. Namely, 
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- statistical physics effects at temperature T, 

- radiative corrections of magnitude M, the spontaneous supersymmetry breaking scale, 

- the presence of moduli fields, whose kinetic energies scale as 

We determine the low energy effective action at the one-loop level and parametrize data at 
the exit Ie oi the Hagedorn phase by a set of IBC. 

To be specific, when n = 1 internal radius is involved in the supersymmetry breaking, we 
show that, depending on the R-symmetry charges of the massless spectra, the models can 
be grouped in different Cases (a), (b) or (c). For each, there exists a partition of the set of 
IBC, whose classes yield cosmological evolutions that share a common behavior. 

• In Case (a), one class of IBC seems to induce a spontaneous decompactification of the 
internal direction involved in the breaking of supersymmetry. It is out of the scope of 
the present paper. 

• In Case (b), the evolutions of the scale factor are monotonic. The expanding solutions 
belong to the same class and are attracted to a radiation dominated era. Asymptoti- 
cally, one has 3H^ oc oc oc M*, where a(t) oc Vt. 

• In Case (c), a class of IBC reaches solutions attracted to a Big Crunch era dominated 
by the radiative corrections to the initially vanishing vacuum energy. Asymptotically, 
one has 3H'^ oc oc l/a^^ oc T^^ where a{t) oc (tec - 1)^^^^- 

• In all Cases, a class of IBC reaches solutions attracted to a Big Crunch era dominated 
by the classical moduli kinetic energy. Asymptotically, one has 3H^ oc oc oc 
]\^i/(i-3/(8o-_))^ where cr_ is given in Eq. (|2.57[) and a(t) oc (tec - tY^^- 



Our analysis has some limitations. First of all, we have supposed that the radii involved 
in the spontaneous breaking of supersymmetry and the radius of the Euclidean time are large 
enough to avoid Hagedorn like singularities in the partition functions. Since we found Big 
Crunch behaviors where T -C M are formally diverging, these evolutions must be restricted 
to M lower than a maximum Hagedorn bound. It would be interesting to relax this constraint 
by considering models free of such singularities [8]. 

Also, we have studied the thermal and quantum corrections at the one-loop level in flat 
backgrounds where moduli have some kinetic energy. This means that the induced sources 
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P, p and their backreactions on the "cold" classical backgrounds we start with must be 
small. Thus, we need afterwords to check the consistency of our approach. 

- For the RDS, we have 

P(^)cxp(t)oc^(x^^O, (5.1) 

which shows its validity. 

- For the Big Crunch attractor a(t) oc (tec — ^)^^^^ of the MDS, both the thermal/quantum 
sources and the classical moduli kinetic energy are diverging. Formally, one has 

= -pW oc oc « » ;^oc^, (6.2) 

which shows that this solution has to be restricted to times where P and p do not exceed 
too much the classical kinetic energy. 

- This limitation docs not exist for the Big Crunch behavior a{t) oc (tec — t)^^'^. Even if 
the thermal/quantum sources and the classical moduli kinetic energy are diverging, one has 

P{t) = -p{t) oc oc ^^^^ oc ^^^^ _ J^3_3/^_ « ^oc^^-^, (5.3) 

since (T_ < 1/2. Thus, the backreaction on the classical background is always small. The 
dynamical attraction to this solution is only limited by the Hagedorn transitions mentioned 
previously. 

To go beyond the constraint that the thermal/quantum corrections must be small, one 
can try to compute the higher string loop corrections and, if possible, use string-string 
dualities to find the non-perturbativc corrections to the effective action. For the models 
under consideration in this work, the extended supersymmetry is only broken spontaneously, 
at the scale M and/or by thermal effects at the scale T. This means that the UV behavior of 
the theories is not affected by these soft breaking and that the renormalization theorems for 
extended supersymmetries should persist. To find the cosmological evolutions beyond the 
first order in perturbation theory, both the thermal effective potential and the corrections 
to the kinetic terms have to be computed. 

The attraction phenomena we found can be generalized to models where n internal radii 
are involved in the spontaneous breaking of supersymmetry. In particular, we have shown 
that for n — 2, the attraction to the radiation era persists. In the process, the complex 
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structure modulus associated to the ratio of the two radii converges to a constant i.e. is 
dynamically stabilized. 

During the convergence to the radiation dominated era, the ratio M/T reaches its critical 
value exponentially or with damping oscillations. Periods of accelerated expansions are also 
allowed. However, the e-fold number seems to depend weakly on n and, at least for n = 1, 2, 
does not give rise to enough periods of inflation (e-fold < 0.2). Since we have considered 
models with A/" = 8,4, 2 initial supersymmetry {i.e. before spontaneous breaking), it would 
be interesting to extend our work to A/" = 1 models and see if the order of magnitude of the 
e-fold is higher. 
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